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13E  TRANSONIC  BUFFET  OF  A SUPERCRITICAL  AIRFOIL 

Reported  by 
X.Z.  Huang 
of  work  by 

B.H.K.  Lee  and  F.C.  Tang,  et  al 


INTRODUCTION 

This  investigation  was  carried  out  in  the  Institute  for  Aerospace  Research  (IAR)  2D  High  Reynolds  Test 
Facility  (Ref.  1 to  Ref.  3 and  Fig.  1)  to  study  the  buffet  characteristics  of  a supercritical  airfoil,  BGK  No.  1 (Fig.  2). 
Steady,  unsteady  surface  pressure  and  normal  force  were  measured  at  various  angles  of  attack  and  Mach  numbers. 
The  statistical  properties  of  the  normal  force  and  pressure  were  carried  out  by  spectral  analyses.  Buffet  onset 
boundaries  were  evaluated  from  the  divergence  of  the  fluctuating  normal  force  while  buffet  intensities  were 
determined  from  the  normal  force  measurements.  The  attached  and  separated  flow  regions  on  the  airfoil  as  well  as 
the  merging  of  a shock  induced  separation  bubble  with  the  trailing  edge  separation  region  were  determined  by  skin 
friction  measurements. 

The  test  program  is  presented  in  Table  1.  There  are  two  BGK  No.l  models.  One  has  normal  static  pressure 
orifices  and  6 pressure  ports  to  measure  pressure  fluctuations  (BGK-1).  Another  has  15  fast  response  transducers 
(BGK-l(m)).  The  model’s  coordinates  and  the  locations  of  pressure  orifices  and  transducers  are  listed  in  Table  2 
Table  3 respectively  (in  CD  ROM).  The  experimental  arrangement  and  results  have  been  described  in  detail  in  Ref. 
4 to  Ref.  9.  Tabulated  data  and  illustrations  are  presented  in  Table  4 lo  Table  7 and  Fig.  3 to  Fig.  16  in  CD  ROM 
with  part  of  the  illustrations  shown  here. 

Fig.  3 and  Table  4 show  the  fluctuating  normal  force  on  BGK-1  model  for  various  Mach  numbers.  Typical 
power  spectra  of  the  normal  force  are  shown  in  Fig.  4.  The  frequencies  of  the  shock  motion  vary  from  70-80  Hz  for 
the  Mach  number  range  of  0.688-0.796  and  arc  partly  listed  in  Table  5.  The  flow  conditions  where  discrete  shock 
oscillations  were  detected  are  summarized  in  Fig.  5.  The  test  program  for  BGK-l(m)  in  Table  1 can  be  sorted  in 
three  cases  as  seen  in  Fig.  5:  1)  points  A,  B,  C,  D,  and  E;  2)  points  a,  b,  c,  d and  e;  and  3)  points  1,  2,  3,  4,  and  5 
respectively.  The  shaded  region  was  obtained  by  fixing  a Mach  number  but  varying  the  incidence  in  the 
experiment.  A power  spectra  plot  of  the  normal  force  was  computed  at  each  a and  the  presence  of  shock  waves 
was  determined  from  observing  whether  the  70-80  Hz  peak  was  present  or  not.  The  buffet  boundary,  which  was 
obtained  from  divergence  of  the  fluctuating  normal  force,  is  included  in  this  figure  for  reference.  This  buffet  onset 
is  identified  from  the  divergence  of  the  normal  force  fluctuations  by  noting  the  point  on  the  curve  with  a slope 
dCrs|/dCL=0.1.  This  value  is  arbitrarily  chosen,  but  in  those  cases  where  buffet  onset  is  primarily  due  to  trailing 
edge  separations,  this  criterion  for  deriving  the  buffet  boundary  is  found  to  give  consistent  results  and  agrees  with 
values  computed  front  trailing  edge  pressure  divergence. 

The  static  surface  pressure  distributions  are  listed  in  Table  6 with  some  examples  shown  here  from  Fig.  6 to 
Fig.  8.  The  cross-hatched  and  open  bar  symbols  in  Fig.  7 and  Fig.  8 denote  regions  of  attached  and  separated  flows 
determined  from  skin  friction  measurements. 

Table  7 presents  the  unsteady  pressure  or  the  pressure  intensities  along  airfoil  chord  of  BGK-1  and  BGK-l(m) 
models.  The  corresponding  figures  are  shown  in  Fig.  10  and  Fig.  1 1 . 

The  statistical  properties  such  as  power  and  cross  power  spectral  density,  auto  and  cross  correlation  functions, 
as  well  as  coherence  functions  of  pressure  and  normal  force  have  been  measured  at  different  Mach  numbers  and 
angles  of  attack.  As  examples  Fig.  12  shows  a set  of  the  spectral  analyses  at  the  condition  of  M=0.753  and  a=5.66° 
for  BGK-1  model.  The  frequency  response  of  the  installed  transducers  was  calibrated  and  established  to  be  flat  up 
to  approximately  200  Hz.  The  normal  force  signal  was  obtained  at  the  sampling  frequency  of  1.6  kHz.  Power 
spectra  of  unsteady  pressure  on  upper  surface  of  BGK-l(m)  at  different  locations  are  shown  from  Fig.  13a  to  Fig. 
13c.  Fig.  14  shows  the  cross  correlation  functions  between  different  transducers  at  M=0.688  and  a=3.99°,  6.43° 
and  9°. 

The  pressure-time  histories  on  BGK-l(m)  model  at  M=0.71  and  various  a are  presented  in  Fig.  15.  The 
unsteady  pressure  fluctuations  behind  the  periodic  shock  wave  have  two  contributions.  One  is  from  a random 
component  associated  with  the  turbulent  motion  in  the  separated  flow  region.  Another  is  a deterministic  part  from 
the  pressure  field  as  a result  of  shock  wave  oscillation.  Thus,  approximately  175  ensemble  averages  of  the  pressure 
signals  were  performed.  Each  ensemble,  which  was  synchronized  to  the  zero  crossings  decided  from  balance 
normal  force  spectra,  had  32  samples.  A Fourier  analysis  was  then  performed  to  obtain  the  fundamental  and 
harmonics  of  the  oscillatory  pressure  field. 
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For  supercritical  airfoils  such  as  the  BGK  No.  1,  it  is  found  that  at  the  lower  Mach  number  range,  separation 
can  occur  behind  the  shock  wave  as  a bubble  and  propagates  downstream  as  the  angle  of  incidence  is  increased. 
Trailing  edge  separation  can  occur  at  the  same  time  and  it  moves  upstream  and  the  two  separated  regions  will 
eventually  merge.  An  investigation  on  the  model  was  carried  out  at  M=0.688  using  a Preston  tube  to  measure  the 
skin  friction  on  the  surface  at  various  angles  of  attack.  The  typical  distributions  of  the  skin  friction  coefficient  are 
presented  in  Fig.  16.  The  results  show  that  at  oc=4.67°,  a small  separation  bubble  begins  to  form  behind  the  shock 
w'ave.  The  separation  bubble  grows  as  the  incidence  increased  and  at  cc=6.15°,  trailing  edge  separation  has  already 
begun  and  has  moved  to  nearly  90%  of  the  chord  as  seen  in  Fig.  16. 


LIST  OF  SYMBOLS  AND  DEFINITIONS 

b model  span 

c model  chord 

CL  lift  coefficient 

Cues  design  lift  coefficient 

Cn  normal  force  coefficient 

Cp  pressure  coefficient 

Cp  ensemble-averaged  pressure  coefficient 

Cp  fluctuating  pressure  coefficient 

C N fluctuating  normal  force  coefficient 

f frequency 

L lift 

M free  stream  Mach  number 

Mdes  design  Mach  number 

Mdr  drag  rise  Mach  number 

N Normal  force 

N time-averaged  normal  force 

Nrms  rms  value  of  normal  force 

P local  static  pressure 

PM  free  stream  static  pressure 

P time-averaged  pressure 

Prms  rms  value  of  the  fluctuating  pressure 

Q.q 


L 

qbc 


N 

qbc 

P~P, 

q 


qbc 


Nrms= 


P = 

1 rms — 


free  stream  dynamic  pressure 
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Re 

Reynolds  number  based  on  chord 

Rx(t) 

auto  correlation  function  of  x(t) 

1 T 

R,(i)  = lim  — Jx(t)  - x(t  + t)dt 
T-»-  T 0 

Rxy(t) 

cross  correlation  function  of  x(t)  and  y(t) 

1 T 

Rxy(t)  = lim  — Jx(t)-  y(t  + T)dt 
T->~  1 0 

Sx(f) 

power  spectral  density  of  x(t) 

Sx(f ) = 2 'jR  X(x)e-i2,tfxdx 

Sxy(f) 

cross  power  spectral  density  of  x(t)  and  y(t) 

Sxy(f)  = 2 JRxy(x)e-i2,I,tdT 

t,T 

time 

X 

distance  measured  along  chord  from  the  leading  edge 

x(t) 

random  signal 

y(t) 

random  signal 

a 

mean  wing  incidence 

Yxy 

coherence  function  of  x(t)  and  y(t) 

.2  |Sxy<r>| 

y Sx(f)  Sy(f) 

1 

time  delay 

FORMULARY 

1 General  Description  of  model 


1.1 

Designation 

Bauer-Garabedian-Korn  (BGK  No.  1)  airfoil 

1.2 

Type 

Aft-loaded,  natural  laminar  flow-capable,  shock-free 
supercritical  airfoil 

1.3 

Design  condition 

Potential  flow 

Mdes=0.72,  MdI=0.75,  ^=0.63 

1.4 

Additional  remarks 

1.5 

References 

Ref.  10 

Model  Geometry 

2.1 

Chord  length 

10  in 

2.2 

Span 

15  in 

2.3 

Model  coordinate 

See  Table  2 in  CD  ROM 

2.4 

Nose  radius 

- 

2.5 

Maximum  thickness 

t/c  = 1 1.8% 

2.6 

Trailing  edge  thickness 

0. 1 % of  the  chord 

2.7 

Additional  remarks 

2.8 

References 

Ref.  4,  10 

Wind  Tunnel 

3.1 

Designation 

IAR  2D  High  Reynolds  Test  Facility 

3.2 

Type  of  tunnel 

Blowdown,  closed  test  section 

3.3 

Test  section  dimensions 

Rectangular,  height  60  in,  width  15  in,  (see  Fig.  la) 
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3.4 

Length  of  parallel  section 

141  in. 

3.5 

Floor  and  ceiling  porosity 

20.5% 

3.6 

Side  wall  boundary  layer 

A gap  between  inlet  and  nozzle  section  permit  bleeding  into 
the  plenum  chamber  of  fairly  thick  side  wall  boundary  layer 
(~2  in.),  see  Fig.  lb. 

3.7 

Side  wall  near  model  area 

Additional  porous  with  boundary  layer  suction  to 
atmospheric,  see  Fig.  lc. 

3.8 

Ventilation  geometry 

See  Fig.  Id. 

3.9 

Range  of  Mach  numbers 

0.1  to  1.1 

3.10 

Re 

40xl06/ft  at  M=l,  10  seconds  total  run  time 

3.11 

Wake  traverse  probe 

7 wafer  (12  ports)  Statham  miniature  transducer  unit 

3.12 

Turbulence  intensity  level 

0.1%  for  Re/ft  < 6x  1 06 
0.16-0.24%  for  Re/ft  10xl06~27xl06 

3.13 

Turbulence  intensity  level 

0. 1%  for  Re/ft  < 6x1 06 
0.16-0.24%  for  Re/ft  10xl06~27xl06 

3.14 

Reference  on  tunnel 

Ref.  1,2  and  3 

Measurements  and  Observations 

4.1 

Steady  pressure  for  the  mean  conditions 

measured  directly 

4.2 

Unsteady  pressure  for  the  mean 
conditions 

measured  directly 

4.3 

Steady  forces  for  the  mean  conditions 

measured  directly 

4.4 

Unsteady  forces  for  the  mean  conditions 

measured  directly 

4.5 

Spectral  analysis  of  the  pressure 

yes 

4.6 

Spectral  analysis  of  the  loads 

yes 

4.7 

Local  skin  friction 

yes 

4.8 

Buffet  boundaries 

yes 

4.9 

Synchronous  Cp  time  histories 

yes 

Test  Conditions 

5.1 

Tunnel  height/model  chord  ratio 

6 

5.2 

Tunnel  width/model  chord  ratio 

1.5 

5.3 

Range  of  Mach  number 

0.501  -0.805 

5.4 

Incidence  range 

-0.36-11.74 

5.5 

Reynolds  number  range 

15xl06  ~20xl05 

5.6 

Range  of  tunnel  total  pressure 

300  psi 

5.7 

Maximum  mass  flow 

10  lbm/sec 

5.8 

Definition  of  model  incidence 

between  “x”  of  model  axis  (Fig.  2)  and  tunnel  axis 

5.9 

Position  of  transition,  if  free 

Not  applicable 

5.10 

Flow  instabilities  during  tests 

No  evidence 

5.11 

Model  deformation  under  the  loads 

Negligible 

5.12 

References  describing  tests 

Ref.  4 to  Ref.  9 

6 Instrumentation 


6,1  Steady  pressure  measurements  for  BGK- 
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1 model 

6.1.1  Position  of  orifices 

6.1.2  Type  of  measuring  system 

6.2  Unsteady  pressure  measurements  for 
BGK-1  model 

6.2. 1 Location  of  transducers 

6.2.2  Type  of  transducers 

6.2.3  Dynamic  response 

6.2.4  Signal  record 

6.2.5  Data  reduction 

6.3  Unsteady  pressure  measurement  for 
model  BGK-1  (m)  model 

6.3.1  Location  of  transducers 

6.3.2  Type  of  transducers 

6.3.3  Diameter  of  screen 

6.3.4  Type  of  screen 

6.3.5  Signal  measurements 

6.3.6  Sampling  rate 

6.4  Loads  measurement 

6.4.1  Type  of  sensors 

6.4.2  Balance 

6.4.3  Pitch  drive  system 


6.4.4  Sampling  rate 
6.5  Skin  friction  measurement 

6.5.1  Type  of  transducers 

6.5.2  Method  of  measurement 

6.5.3  Spatial  resolution 

7 Data  presentation 

7.1  Test  cases 

7.2  Normal  force  fluctuation 

7.3  Shock  oscillation  frequencies 

7.4  Region  of  shock  oscillation 

7.5  Steady  pressure 

7.6  Unsteady  pressure 

7.7  Spectral  analysis 

7.7.1  Power  spectral  density 


See  Fig.  2a  and  Table  3 in  CD  ROM 
70  pressure  tubes  + 1 5 in  situ  pressure  transducers 

See  Fig.  2a  (in  the  middle  chord)  and  Table  3 in  CD  ROM 
6 Kulite  TQ  360  25  psid  transducers 
Flat  up  to  approximately  200  Hz. 

Recorded  on  FM  tape  for  subsequent  analysis. 


See  Fig.  2b  and  Table  3 in  CD  ROM 

16  of  25  psid  custom  made  CQ-062-25D  differential  Kulite 
transducers 

0.042  in. 

0.005  in  thick  with  0.062  in  diameter  holes  in  a mesh 
pattern 

Signals  were  filtered  by  a four  pole  low  pass  filter  having  a 
300  Hz  3db  point  and  a -24  db/octave  slope  beyond  600 
Hz. 

1.6  kHz 
strain  gages 

3 component  side  balance  with  max  capacity  of  N-20.000 
lbf,  m= 22,500  in.lb  and  X=2,000  lbf 

Range:  55° 

maximum  angular  rate:  12°  /sec,  fully  loaded 
step  program:  0.25°,  0.5°,  1°,  2°,  5° 
ramp  program:  0°  - 10°/sec 
1.6  kHz 

Given  by  the  difference  between  the  total  and  static 
pressures 

Preston  tube  to  determine  the  pitot  pressure 
0.05c  for  x > 0.6c  and  0.02c  for  x < 0.6c  respectively 

See  Table  1 

Fig.  3,  Fig.  4 and  Table  4 in  CD  ROM 
Table  5 in  CD  ROM 
Fig.  5 

Fig.  6 to  Fig.  9 and  Table  6 in  CD  ROM 
Fig.  10,  Fig.  1 1 and  Table  7 in  CD  ROM 
Fig.  12,  Fig.  13  and  Fig.  14 
Fig.  12a  and  Fig.  13 


324 


7.7.2  Auto  correlation  functions 

Fig.  12b 

7.7.3  Cross  correlation  functions 

Fig.  12c  and  Fig.  14 

7.7.4  Cross  power  spectral  density 

Fig.  1 2d 

7.7.5  Coherence  function 

Fig.  12e 

7.7.6  Cross  power  spectral  density  and 

Fig.  12f 

coherence  function  between  pressure 
and  normal  force 

7.7.7  Pressure-time  histories 

Fig.  15 

7.8  Skin  friction 

Fig.  16 

7.9  Example  illustrations  of  results 

Fig.  6,  Fig.  7,  Fig.  8,  Fig.  10  to  Fig.  16 

Comments  on  data 

8.1  Mach  number 

Mach  number  be  maintained  constant  by  conLrol  system 

8.2  Steady  incidence 

measured  by  a potentiometer 

8.3  Balance  linearity 

maximum  0.3%  and  generally  <0.1% 

8.4  Balance  interaction 

<1.26% 

8.5  Balance  natural  frequencies 

140,  215,  320,  360  Hz,  buffet  excitation  frequencies=70-80 
Hz 

8.6  Unsteady  pressure  coefficients 

a discrete  frequency  of  *420  Hz  was  detected  due  to  tunnel 
disturbances  (See  Fig.  4) 

8.7  Wall  interference  corrections 

distributed  suction  was  applied  through  porous  plates  in  the 

vicinity  of  the  model  to  minimize  any  three-dimensional 
effects 


9 Personal  contact  for  further  information 

X.Z.  Huang,  Aerodynamics  Laboratory,  Institute  for  Aerospace  Research,  National  Research  Council  of 
Canada 

M-10,  IAR/NRC,  Montreal  Rd.  Ottawa,  Ontario,  Canada,  K1A  0R6 
e-mail:  xingzhong.huang@nrc.ca 
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Pressure 


2b  BGK-l(m)  model 

Fig.  2 BGK  No.  1 supercritical  airfoil 


Fig.  5 Region  of  shock  oscillation 
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Fig.  12d  (cont)  Cross  power 
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Fig.  12e  Coherence  functions  for  M„=0.753,  CL=0.945,  q=24.5  psi 
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Fig.  12f  Cross  power  spectral  density  and  coherence  function  between  pressure  and 
normal  force  for  M„=0.753,  CL=0.945,  q=24.5  psi 
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